Reaction-diffusion processes of three species on small-world networks 
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We study the decay process for the reaction-diffusion process of three species on the small-world 
network. The decay process is manipulated from the deterministic rate equation of three species in 
the reaction-diffusion system. The particle density and the global reaction rate on a two dimensional 
small-world network adding new random links is discussed numerically, and the global reaction 
rate before and after the crossover is also found by means of the Monte Carlo simulation. The 
time- dependent global reaction rate scales as a power law with the scaling exponent 0.66 at early 
time regime while it scales with —0.50 at long time regime, in all four cases of the added probability 
p = 0.2 — 0.8. Especially, our result presented is compared with the numerical calculation of regular 
networks. 



I. INTRODUCTION 

In the past decades, several papers have been devoted 
to the theoretical and numerical investigation on the re- 
actant segregation phenomena in anomalous kinetics. 1-3 
An important contribution to this field of research has 
been given by Ovchinnikov and Zeldovich, 4 who devel- 
oped the segregation phenomenon. It is well known that 
the binary reaction has primarily been investigated on 
the process of ternary reactions 5 and many phenomena 
in nature. 6-10 Until now, the segregation of reactants 
has been investigated a reaction-diffusion process based 
on assumption that the scaling form for A + B — > C 
has almost been verified by the computer simulation and 
the experiment. 11-15 The decay process for the reaction- 
diffusion system have been investigated by Taitelbaum et 
al. 16 who have precisely shown both the global reaction 
rate and the reaction front increased as a function of i 1 / 2 
at very early times. At that time, Cornell et al. 17,18 have 
argued the diffusion-limited reaction nA + mB — > C for 
both homogeneous and inhomogeneous conditions under 
initially separated reactants. They have discussed that 
the global reaction rate decreases as t -1 ' 2 in long time 
limit, independent of n and m, and that the upper critical 
dimension is d c = 2 for the reaction-diffusion process. 

Recently, Zumofen et al. 19 have studied mainly the 
particle density and the pair correlation function on the 
two-particles reaction process. Particularly, they found 
that the particle density is distributed as P(n) ~ n -1-7 , 
n > 0, 1 < 7 < 2, for Levy walks and that the segregation 
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disappears in d = 3 dimension for 7 < 3/2. Yen et al. 20 
have studied for the asymptotic early-time scaling in the 
ternary reaction-diffusion process with initially separated 
reactants. Moreover, Kim et al. 21 argued the decay pro- 
cess in a reaction-diffusion system with three species on 
the regular square lattice. By a simple perturbation ex- 
pansion, we have analytically derived the particle density 
and the global reaction rate before and after the crossover 
in the reaction-diffusion system of A + B + C — > 0. 

For last few years, a growing interest has been concen- 
trated on small-world and scale-free network models, 22 
and these models have recently been widely studied in 
various applications of methods from physics to other 
natural, social and applied sciences. In fact, the phenom- 
ena of small-world and scale-free networks are showed to 
be different from the dynamical behavior of the regular 
lattice system. The degree distribution for scale-free net- 
works scales as a power law p(k) ~ fc -7 for large k, while 
it decays faster than exponentially for random networks. 
It is, in present, of fundamental importance that the nu- 
merical and analytical result of small-world and scale- free 
network models is compared with that of regular network 
models. 

Recently, the bimolecular chemical reaction in scale- 
free networks 23 is in scale-free networks studied the gen- 
eration of the depletion zone and the segregation of the 
reactants. It was found that the reaction-diffusion pro- 
cesses in scale-free networks is different in their nature 
compared to regular lattice models, due to the small di- 
ameter of networks and the existence of hubs. Further- 
more, Catanzaro et al. have analyzed that the inverse 
particle density scales linearly as 1/pit) ~ t. From this 
result, they have found that the inverse particle density 
in uncorrelated scale- free network crosses over to a linear 
behavior. Very recently, Gallos and Argyrakis 25 have dis- 
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cussed the reaction-diffusion process of the two species on 
the scale-free network between the the correlated and the 
uncorrelated configuration models. They have especially 
revealed that two models are identical when 7 = 3.0. 

In this paper, the decay process for the reaction- 
diffusion process of three species on the small-world net- 
work is studied. We also consider the particle density and 
the global reaction rate on the two-dimensional small- 
world lattices added by new random links. In Section 2, 
we discuss the reaction-diffusion process of three species 
on small- world networks. We present some results ob- 
tained by the numerical simulations and the concluding 
remark in the final section. 



II. REACTION-DIFFUSION SYSTEM OF 
THREE SPECIES 

In the reactions of A + A — > and A + B — > types 
on a d dimensional regular lattice, the surviving particle 
density scales in long time limit as 

4t- — ~* q , (1) 
p(t) Po V ^ 

where p is the initial particle density. It is well known 
that the scaling exponent a = d/d c for d < d c and 1 for 
d> d c , and that the critical dimension d c = 2 for A + A 
and 4 for A + B. 

Next the A + B — > process on a regular network^ 19) 
is introduced for Levy walks using the following reaction- 
diffusion equation, 

JU(f, t) = DLA(r, t) - nA(r, t)B(r, t), (2) 

-B(r, t) = DLB{r, t) - nA(f, t)B{r, t), (3) 

where A(f, t) and B(r,t) are the particle densities, D a 
generalized diffusion coefficient, 7 a reaction rate, and L 
the operator for the Levy-enhanced diffusion. The time- 
dependent particle densities can be calculated as 

A{t) = B{t) ~ t d ' 2 -<, for 7 > d/2, (4) 

where three marginal values are 7 = 1 for d = 2, 7 = 3/2 
for d — 3, and 7 = 2 for d = 4. 

From the decay process 20 of A + 2B^C, the diffusion 
equations under the initial reactant segregation are as 
follows: 

%-A(x,t) = D a V 2 A(x,t) -kA(x,t)B 2 (x,t) 1 (5) 



— B(x,t) = D b X7 2 B(x,t) - kA(x,t)B 2 (x,t), (6) 

where A(x,t) and B(x,t) are the particle densities, k is 
the microscopic reaction constant, and D a and Df, are 



the diffusion coefficients of reactants A(x,t) and B(x,t), 
respectively. Then, from the lowest order of the pertur- 
bation theory, the global reaction rates R(t) on both the 
early time and long time behaviors scaled as a power 
law t 1 / 2 and i -1 / 2 , respectively. Using the Monte carlo 
method, it was found that the slopes are, respectively, 
0.5 and —0.48 before and after the crossover. 

Let us denote that A(x, t), B(x, t), and C(x, t) are the 
particle densities for three-species A, B, and C existing at 
a position x at time t. We assume that three species are 
initially distributed separately on the axis x. Then, the 
deterministic rate equation for A(x, t) can be expressed 
in terms of 

— A(x, t) = D A V 2 A(x, t) - KA(x, t)B(x, t)C(x, t), (7) 

where Da is the diffusion constant for one species A, and 
K is the reaction rate. The solution for A(x,t), B(x,t), 
and C(x,t) is obtained that 

A(x, t) - B(x, t) - C(x, t) ~ 0( J!L) (8) 

in large time limit, where (f>(x) = ^= dxe~ x . Sim- 
ilar to the result of the previous work, 20 the time de- 
pendence of the global reaction rates R(t) was shown to 
behave as 

R(t) - t 1 ' 2 and r 1 ' 2 (9) 
in early time and long time limits. 21 



III. NUMERICAL RESULTS AND 
CONCLUDING REMARKS 

In order to confirm numerically the global reaction 
rate, we mainly use the Monte Carlo simulation method 
in reaction-diffusion process of^l + B + C^Oon the 
small- world network adding new random links. First of 
all, we assume that three species reactants are distributed 
randomly on two dimensional square lattice with the pe- 
riodic boundary condition. After one species reactant 
is chosen at random, the direction of its move is cho- 
sen at random with equal probability to one of its linked 
neighbor lattice points. When two species reactants meet 
each other on the same lattice point, the intermediate 
process existing concurrently the combined two species 
reactants can be formed. If two species reactants meet 
the third reactant, these reactants react and leave imme- 
diately on lattice point. If one species or two species re- 
actants cannot meet the other reactant, the reactants dif- 
fuse randomly to one of its linked neighbor lattice points. 
The diffusion constant in our case takes the same value 
for each reactant, and the respective particle density of 
13.33(40/3)% for A, B, and C is distributed randomly on 
a square lattice having 200 x 200 lattice points with the 
boundary condition. After our simulation is performed 
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FIG. 1: Plot of the global reaction rate R(t) for the reaction- 
diffusion process of A + B + C — > 0. Our numerical simulation 
is performed on 2 x 10 3 configurations in 200 x 200 square 
lattice that is added by new random links for the added prob- 
ability p — (black circles), 0.2 (squares), 0.4 (circles), 0.6 
(down triangles), and 0.8 (up triangles). 

on 2 x 10 3 realizations and K = 1/1500, we directly ob- 
served the crossover for the global reaction rate from our 
simulation result. It is numerically found that the scal- 



ing exponent for the slope before and after the crossover 
is, respectively, 0.66 and —0.54 on the regular network 
(p = 0). As shown in Fig. 1, the time-dependent global 
reaction rate scales as a power law with the scaling expo- 
nent 0.66 at early time regime while it scales with —0.50 
at long time regime, in all four cases of the added prob- 
ability p = 0.2, 0.4, 0.6, and 0.8. 

In conclusion, we have numerically estimated the 
global reaction rate before and after the crossover in 
reaction-diffusion system of A + B + C — » on square 
lattice added by new random links. It is really found 
from our simulation result that the scaling exponent of 
the global reaction rate on small-world network added 
by new random links is the similar to that of the regular 
network at early time regime. At long time regime, the 
decay process on small-world network proceeds slightly 
faster than in the case of the regular network. 

In future, our work is in progress to extend the cor- 
related and the uncorrelated configuration models. We 
also will attempt to investigate small-world and scale- 
free networks in several scientific fields for three species 
reaction-diffusion process. 
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